Studies of Thermal Conductivity in Carbon Nanotubes Using Simulation
Sunil H. Ganatra*, Sneha D. Khobragade, Anushree S. Ujjankar, Chudaman D. Pourkar
Department of Chemistry, Institute of Science, R. T. Road, Nagpur – 440001, M.S. India
*Corresponding Author E-mail: sunilganatra@gmail.com
ABSTRACT:
The Nano Heat flow tool allows users to explore the time evolution of kinetic and potential energy among the vibrational modes of a carbon nanotube over the course of a molecular dynamics (MD) simulation. It is possible to observe the drop of vibrational energy through the modes of the system as a non-equilibrium population of phonons is dissipated towards thermal equilibrium, and thus gives insight into the basic sources of damping and dissipation within nanoscale objects.
Molecular Dynamics Simulation for Heat Flow in C(6,0), C(6,3) and C(6,6) carbon nanotubes were studied at 289.15 K using high performance parallel processing computers.
It is reported that the selected nanotubes i.e. C(6,0), C(6,3) and C(6,6) shows the thermal conductivity at various frequencies.
The surprise results for C(6,0) and C(6,6) where the frequency is highly selective and give single peak. Whereas, in case of C(6,3) various frequency peaks are reported and spread in wide spectrum.
KEYWORDS: Carbon nanotubes, Heat Flow, Thermal conductivity, Molecular dynamics.
Carbon nanotubes were discovered in 1991 by Iijima of NEC Corporation.[1] Since then, efforts in synthesis, characterization and theoretical investigation on nanotubes has grown exponentially. This is due to their novel mechanical, electronic properties and tremendous potential for future technological applications. In 1993, the simplest kind of carbon nanotube i.e. single walled carbon nanotubes (SWNTs) were discovered independently by Iijima group and an IBM team headed by Bethune [2]. These SWNTs can be regarded as a rolled-up graphite sheet in the cylindrical form. Some specific defect-free forms of these SWNT show remarkable mechanical properties and metallic behaviour [1,2].
The Nano Heat Flow Tool [3] allows users to explore the time evolution of kinetic and potential energy among the vibrational modes of a carbon nanotube over the course of a molecular dynamics (MD) simulation. It is possible to observe the drop of vibrational energy through the modes of the system as a non-equilibrium population of phonons is dissipated towards thermal equilibrium, and thus gives insight into the basic sources of damping and dissipation within nanoscale objects.
Figure 1: Nano Vibrational Modes in Swnts.
Material and Methods:
Molecular Dynamics calculations were performed using Large-scale atomic/molecular massively parallel simulator LAMMPS on carbon nanotubes using online parallel processing computers (nanohub.org) [3] to understand the heat flow in C(6,0), C(6,3) and C(6,6). The nanotube coordinates were generated by Tube Gen Online software [4,5]. Micro canonical ensemble used to perform classical MD simulation of selected nanotubes. The initial conditions for the experiment restricted to the excitation of single vibrational modes of low frequency [6,7,8].
The simulations were performed on NanoHub [3] with Nano Heat Flow Tools. All tubes were tested for heat flow along the tube using LAMMPS. Simulations for all three tubes were performed separately with same basic parameters. [9,10,11] as listed in Table 1.
Table 1. : MD simulation Input Parameters
Parameter Name |
Value |
Remark |
|
Type of Nanotube |
C(6,0) : Zigzag |
One type at a time. In total three simulations performed for each type of tube. |
|
Nanotube Unit Cells |
2 |
|
|
Temperature |
298.15 K |
|
|
Number of Steps |
20000 |
Number of steps of MD simulation. |
|
Step Size |
0.2 fs |
Time period for each step size. |
|
Vibrational Mode |
First Longitudinal |
|
Table 2 : The peak of energy / frequency graph for C(6,0) type Carbon Nanotube
|
Sr.No. |
Energy |
Frequency |
|
1 |
0.060919 |
22.39164 |
|
2 |
0.080844 |
22.72093 |
|
3 |
0.11228 |
23.05022 |
|
4 |
0.165961 |
23.37951 |
|
5 |
0.268303 |
23.7088 |
|
6 |
0.497221 |
24.03809 |
|
7 |
1.139402 |
24.36738 |
|
8 |
3.062294 |
24.69667 |
|
9 |
2.794253 |
25.02596 |
|
10 |
1.0292 |
25.35524 |
|
11 |
0.461296 |
25.68453 |
|
12 |
0.253395 |
26.01382 |
|
13 |
0.158543 |
26.34311 |
|
14 |
0.108098 |
26.6724 |
|
15 |
0.078267 |
27.00169 |
|
16 |
0.059224 |
27.33098 |
|
17 |
0.046348 |
27.66027 |
|
18 |
0.037245 |
27.98956 |
Table 3 : The peak of energy / frequency graph for C(6,6) type Carbon Nanotube
|
Sr.No. |
Energy |
Frequency |
|
1 |
0.05045 |
36.27509 |
|
2 |
0.06663058 |
36.63073 |
|
3 |
0.09196963 |
36.98637 |
|
4 |
0.1348236 |
37.34201 |
|
5 |
0.2155027 |
37.69765 |
|
6 |
0.3932584 |
38.05328 |
|
7 |
0.8891345 |
38.40892 |
|
8 |
2.591741 |
38.76456 |
|
9 |
3.064494 |
39.1202 |
|
10 |
1.056946 |
39.47584 |
|
11 |
0.4453835 |
39.83148 |
|
12 |
0.2367583 |
40.18711 |
|
13 |
0.1453171 |
40.54275 |
|
14 |
0.09786072 |
40.89839 |
|
15 |
0.070249 |
41.25403 |
|
16 |
0.052824 |
41.60967 |
|
17 |
0.041143 |
41.9653 |
|
18 |
0.032938 |
42.32094 |
Table 4 : The peak of energy / frequency graph for C(6,3) type Carbon Nanotube
|
Sr.No. |
Energy |
Frequency |
Frequency
Range |
|
1 |
3.606208 |
19.77989 |
16.6567 –22.9030 |
|
2 |
1.992325 |
9.716438 |
7.2872 – 11.1045 |
|
3 |
1.305825 |
39.21277 |
36.7836 – 41.6418 |
|
4 |
1.136508 |
28.45528 |
26.7202 –30.5373 |
|
5 |
0.7431926 |
56.91056 |
55.5225 – 57.6047 |
|
6 |
0.7148732 |
58.29863 |
57.6046 – 60.3807 |
|
7 |
0.403017 |
49.27622 |
47.8881 – 50.6642 |
|
8 |
0.308775 |
12.1455 |
11.4515 – 13.1865 |
|
9 |
0.2553308 |
32.27245 |
31.2314 – 32.9664 |
|
10 |
0.2066284 |
66.627 |
65.5859 – 69.0561 |
|
11 |
0.1551181 |
54.82847 |
53.7874 – 55.5225 |
|
12 |
0.1447582 |
61.07475 |
60.0337 – 62.1158 |
|
13 |
0.1147414 |
64.19789 |
63.5038 – 65.2389 |
|
14 |
0.111506 |
52.39936 |
51.7053 – 53.4401 |
Results:
The MD simulations were executed and the heat flow values in the form of energy (on arbitrary scale) against the frequency were calculated. The reported values for C(6,0), C(6,6) and C(6,3) nanotubes are shown in Table 2,3 and 4 respectively. Table 5 shows the highest energy values at specific frequencies for C(6,0), C(6,6) and C(6,3) types of Carbon nanotubes at 298.15 K.
The variations of Eigen value frequency in THz were studied for various vibration modes for all three types of carbon nanotubes. Figure 3 shows the variation of Eigen value frequency in graphical form. It was clearly observed that C(6,0) and C(6,6) represents the uniform phenomenon, whereas C(6,3) shows highly diverse phenomenon.
Figure 4, 5 and 6 shows the 3D graphical representation of time scale variation of Energy Vs Frequency for C(6,0), C(6,6) and C(6,3) Carbon nanotube.
Table 5 : The Highest Energy Values at a specific frequency for various types of Carbon Nanotubes at 298.15 K.
|
Sr. No. |
Nanotube Type |
No. of Peaks |
Highest Energy Peak |
Frequency (THz) |
Range of Frequency for peak energy |
|
1 |
C(6,0) |
1 |
3.062294 |
24.69667 |
22.39164– 27.98956 |
|
2 |
C(6,6) |
1 |
3.064494 |
39.1202 |
36.27509– 42.32094 |
|
3 |
C(6,3) |
14 |
3.606208 |
19.77989 |
16.6567 –22.9030 |
Figure 2 shows the variation of frequencies with energy for C(6,0), C(6,6) and C(6,3) types of Carbon nanotubes at 298.15 K.
Figure 2 : Variation of frequencies with energy at 298.15 K.
fig. 3 : variation of Eigen value verses vibrational mode number.
fig. 4 : 3D Time scale variation of energy vs frequency for c(6,0)
Fig. 5 : 3D Time Scale Variation of Energy Verses Frequency For C(6,6)
Fig. 6 : 3D Time Scale Variation of Energy Verses Frequency For C(6,3)
Conclusion:
The selected nanotubes i.e. C(6,0), C(6,3) and C(6,6) show the thermal conductivity at various frequencies. The surprise results for C(6,0) and C(6,6); where the frequency is highly selective and having single peak. Whereas, in case of C(6,3) various frequency peaks are reported and spread in wide spectrum.
The C(6,0) and C(6,6) shows heat energy at 25 THz and 38 THz respectively. Only single peak is observed in these types of nanotubes. This uniqueness is due to their structure arrangement. Whereas C(6,3) shows heat energy at various levels. There are number of peaks at 10, 20, 28, 38, 48 THz frequency with single peak and 58 THz with shoulder. This is again due to the structure arrangement of carbon in C(6,3) structure. The heat flow in this case of carbon tube is with multiple frequencies and spread widely on frequency band.
The conclusion is also confirmed from the figure 4,5 and 6 showing the time scale variation of Energy Vs Frequency. In case of C(6,0) and C(6,6) the frequency is uniformed whereas in C(6,3) the frequency is varied with energy.
REFERENCE:
1. Iijima, S., Helical microtubules of graphitic carbon," Nature, 354:1991: 56-58.
2. Iijima, S. and Ichlhashi, T., Single-shell carbon nanotubes of 1-nm diameter, Nature, 363:1993: 603-605.
3. Joe Ringgenberg et al., Nano Heatflow, DOI: 10254/nanohub-r3307.4. , 2007
4. N. Mingo and D. A. Broido, Carbon Nanotube Ballistic Thermal Conductance and Its Limits, DOI: 10.1103/PhysRevLett.95.096105.
5. L. G. C. Rego and G. Kirczenow, Phys. Rev. Lett. 81:1998: 232.
6. M. S. Dresselhaus and P. C. Eklund, Adv. Phys. 49, 200:705.
7. V. P. Carey et al., Review of Heat Transfer Physics, Nanoscale and Microscale Thermophysical Engineering, DOI:10.1080/15567260801917520
8. S. Berber, Y.-K. Kwon, and D. Toma´nek, Phys. Rev. Lett. 84: 2000: 4613.
9. J. Hone, M. Whitney, C. Piskoti, and A. Zettl, Phys. Rev. B 59: 1999:R2514 1999.
10. C. Yu, L. Shi, Z. Yao, D. Li, and A. Majumdar, Nano Lett. 5, 2005:1842 .
11. Bethune, D.S. et.al., "Cobalt-catalyzed growth of carbon nanotubes with single- atomic-layer walls," Nature, 363:1993: 605-607.
Received on 28.03.2012 Modified on 05.04.2012
Accepted on 15.04.2012 © AJRC All right reserved
Asian J. Research Chem. 5(4): April 2012; Page 500-503